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Abstract 

We present thermal conductivity measurements performed in three short-period (GaAs)g(AlAs)5 
superlattices. The samples were grown at different temperatures, leading to different small scale 
roughness and broadening of the interfaces. The cross-plane conductivity is measured with a 
differential 3uj method, at room temperature. The order of magnitude of the overall thermal 
conductivity variation is consistent with existing theoretical models, although the actual variation 
is smaller than expected. 

PACS numbers: 66.70.-f, 68.65.-k, 68.65.Cd, 44.10.+i 
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I. INTRODUCTION 



The thermal conductivity of semiconductor superlattices is strongly reduced with respect 
to the bulk values of their constituents Q, \l , . The heat is mainly carried by phonons 

and two mechanisms can explain the effects of the nanostructuration on thermal transport. 
According to the first (intrinsic) mechanism, the zone- folding involves a modification of the 

n q 

phonon dispersion and allows Umklapp process for low energy phonons [6J, [7fl. According to 
;he second {extrinsic) mechanism, phonons are scattered by imperfections at the interfaces 
8, 9]. Of course, both mechanisms can be active but their relative contributions is still under 
debate, mainly for short period superlattices [loj]. The assessement of the interface role relies 
on discrepancies between theoretical calculations, based on intrinsic mechanisms only, anc 
experimental observations. In order to take into account interfaces defects, some authors 9j 
combined the phonon dispersion curves arising from zone folding and interface scattering in 
order to fit experimental results on GaAs/AlAs superlattices |3|. Mini-Umklapp processes 
were not considered and the phonon mean free path in the layers was supposed to be identi- 
cal to the corresponding bulk value. These authors derived that 17% of the incident phonons 
on an interface undergo diffuse scattering, because of interface defects. Recently, other au- 
thors solved the Boltzmann's equation in superlattices, beyong the constant relaxation time 
approximation [10]. They adressed the role of the mini-Umklapp process and computed 
the contribution of the intrinsic mechanism to the thermal conductivity. Comparing their 
results with experimental data, they inferred that, in (GaAs)s(AlAs)3, "the reduction due 
to extrinsic phonon scattering is roughly three times larger than that due to intrinsic scat- 
tering". Molecular dynamics calculations also stressed the importance of interface defects 
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12J. GaAs/AlAs superlattices were simulated with a simplified structure: the two-atom 



unit cells of GaAs and AlAs are substituted for single average atoms and the interface de- 
fects are obtained by assigning at random the atomic sites in the last monoatomic layer 
of each superlattice layer to an (average) atom or to the other, with a given probability. 
A 60% decrease of the thermal conductivity was then predicted in (GaAs)y(AlAs)7 for a 



50% substitution probability |ll|, ll2j. Lastly, the temperature dependence of the thermal 



conductivity gives clues on interface scattering effects [3|. 

To our knowledge, no experiment has addressed directly the role of interface imperfec- 
tions. In this paper, we report on experiments performed in (GaAs)g(AlAs)5 superlattices 
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exhibiting different interface width. We have measured their cross-plane thermal conductiv- 
ity, using the so-called differential 3uj technique, at room temperature jisl ]. 

II. EXPERIMENTAL RESULTS 

The samples are GaAs/AlAs superlattices grown by molecular beam epitaxy, on GaAs 
buffers. The nominal widths of GaAs and AlAs layers are 2.5 nm and 1.5 nm, respec- 
tively, and the number of GaAs-AlAs periods is 120. Table (111) displays the actual periods 
determined by X-ray diffraction, versus samples. 

Ideal interfaces are infinite atomically flat planes separating pure GaAs and AlAs layers. 
However, actual interfaces exhibit islandlike structures, characterized by their lateral extent 
and by their height (in terms of monolayers). Direct interfaces (AlAs grown on GaAs) tend 
to exhibit large islands (lateral extent from ~ 100 to ~ 1000 nm), whereas inverse inter- 
faces (GaAs grown on AlAs) tend to exhibit smaller islands (lateral extent ~ 0.1 to ~ 10 
nm). Compositional fluctuations inside large islands may also occur. The quality of actual 
interfaces depend on a number of parameters, including growth temperature, interruption 



time... and are usually characterized either by photoluminescence 
tering 
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14| or by Raman scat- 



161 ] . Photoluminescence is sensitive to large scale fluctuations, whereas Raman 



scattering is sensitive to small scale fluctuations. Small scale fluctuations lead to an effective 



concentration profile at the interfaces which can 
be used to characterized the interface width 16 



3e evaluated by Raman scattering and can 



Our samples were grown under the same conditions, except for the substrate temperature 



T s , between 510 and 650 C. 
to the growth temperature 



n that range, small scale interface fluctuations are very sensitive 



15| 



Indeed, the very samples that we have studied were characterized previously by Raman 



scattering [16] . The intensity and frequency shift of the Raman lines were analyzed in term 
of interface broadening. Identical broadening was assumed for all interfaces. Table ([TJ) gives 
the interface width d , determined from the analysis of the confined optical modes. The 



samples are labelled with the same name SI, S2 and S4 as in [16]. Clearly, the interface 
width increases as the growth temperature increases. 

n 

We use the differential 3u method [13J to measure thermal conductivity of three superlat- 
tices, grown on GaAs substrates: SI, S2 and S4. We compare their thermal responses with 
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the responses of GaAs test samples Tl, T2 and T3. Previously, a thin (100 nm) dielectric 
SiC>2 layer is sputtered on SI, S2, S4, Tl, T2 and T3 in order to insure electrical isolation 



between thermal transducers and samples 18]. The insulating SiC>2 layer is deposited in 



a same run on the samples and test samples. Moreover, during sputtering, the samples 
and test samples are mounted on a rotating substrate holder. In this way, the insulating 
layers thickness is found to be nearly constant on Tl, T2 and T3: 104, 102 and 105 nm, 
as measured by ellipsometry. We extrapolate that insulating layer thickness is constant on 
SI, S2, S4, Tl, T2 and T3, within 1.5%. Thermal transducers are processed by lift-off 
photolithography and thermal evaporation of gold (200 nm) on a thin chromium adhesion 
layer (few nm). The typical line width is 30 /im (SI: 35.9 pm, S2: 35.7 pm, S3: 33.1 pm, 
Tl: 28.2 /im, T2: 27.3 /im, T3: 28.0 pm). The line length is 2.5 mm (SI, S2, S4) or 
3.0 mm (Tl, T2, T3). The ratio of the heater line width to the films thickness is around 
50. This ratio is large enough to insure one-dimensional heat flow through the films and 



measurement of cross-plane conductivities. Indeed, simulations based on |19j] confirm that 
possible anisotropy of the superlattices has negligeable effect on the cross-plane conductivity 
measurement. 

The heater /thermometer line thermal response T 2uj , i.e. the in-phase and quadrature 
components of the temperature oscillation at 2u } reads [3] : 

r> /*+oo • 2 

Ti^JL -^L= iy + AT (1) 
IkKJq y 2 ^y 2 + m 2 

where : 

* = 4 < 2 » 



a = (3) 

pCw 2 



AT = (4) 



At low frequency, equation ([T]) can be approximated by |l8| : 
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+ AT (5) 



uj/2n is current frequency supplied to the line. P is the power supplied to the line at 2uj. 
A, p and C are the thermal conductivity, mass density and mass specific heat of the GaAs 
buffer. AT and R c are the temperature drop and thermal resistance through the insulating 




FIG. 1: (Color online) Typical in-phase and out-of-phase temperature oscillations of the 
heater/thermometer line, versus current frequency. Input power: 3 W.m . Symbols: experi- 
mental points. Lines: fits using eq. ([!]). 
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FIG. 2: (Color online) Fitting parameter AT (eq.([T])) versus power flux through the 
heater /thermometer line. Symbols: experimental points. Lines: linear fits. 

layer (Tl, T2, T3) or through the insulating layer and superlattice (SI, S2, S4). I and w are 
the length and width of the heater /thermometer line. K ~ 0.9066 is a constant parameter. 

The thermal response of the samples and test samples is recorded at 21°C, versus fre- 
quency for P/l ranging from 0.5 to 3.0 W m _1 . We fit experimental data (both in-phase and 
quadrature components), using equation ([I]) and two fitting parameters: A and AT. All the 
other terms are either measured (I, w, P) or extracted from litterature (p = 5317 kg m~ 3 , 
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FIG. 3: Thermal conductivity of (GaAs)g(AlAs)5 superlattices, versus interfaces width. 



TABLE I: Sample parameters and experimental results. The superlattices were grown at a substrate 
temperature of T s . Their period, interface width [16] and thickness are d, do and e, respectively. R c 
is the thermal resistance through SiC>2 layers (Tl, T2, T3) or through SiC>2 layers and superlattices 
(SI, S2, S4). R s i and A s ; are the superlattice thermal resistance and thermal conductivity. 



Sample 


T s d 


d 


e 


10 7 x R c 10 7 x R sl A sl 




(°C) (nm) 


(nm) 


(nm) 


(W- x m 2 K) (W~ l m 2 K) (Wm^K' 1 ) 


Tl 








1.18 


T2 








1.16 


T3 








1.21 


SI 


510 4.01 





481.2 


1.74 0.55 8.7 


S2 


550 4.01 


0.15 


481.2 


1.69 0.50 9.5 


S4 


650 3.90 


0.45 


468.0 


1.95 0.76 6.2 



C = 326 J K^kg- 1 20]). The uncertainty on p and C does not affect our conclusions 
because the comparison between samples and test samples eliminates the substrate con- 
tribution. Figure ([I]) displays the experimental thermal response T 2uJ and fits obtained at 
P/l = 3 W m _1 on SI, S2 and S4. Figure is a plot of the fitting parameter AT versus 
the input power flux P/(lw). For a given sample, the points are aligned and, according to 
equation (j3J), the slope is the thermal resistance R c . Table ([!]) summarizes our results. 



6 



The thermal resistance R c through the insulating layer and superlattice (SI, S2, S4) is 
simply: 



R c = Ro + R s i (6) 
Rsi = {- (7) 

where Ro and R s i are the thermal resistances of the insulating SiC>2 layer and of the su- 
perlattice. e and A s ; are the superlattice thickness and thermal conductivity. In order to 
obtain R s i, we use Ro = 1.19 10~ 7 W^ 1 m 2 K, which is the value we measure on Tl, T2, 
T3 within ±2%. Table ((T|) summarizes our results. Our results lie in the same range as 
previously reported for GaAs/AlAs superlattices, although direct comparison is hindered by 
the different periodicities of the samples 3]]. 

The uncertainty on R c comes mainly from uncertainty on width w of the 
heater /thermometer line. We estimate Aw = 0.5/xm. This induces AR C /R C = 3.2%. From 
those figures, we deduce that layer thermal resistances are nearly equal for SI and S2 and 
significantly larger on S4. Actually, we are interested in the variation of A s ; versus samples. 
As long as those variations are concerned, Ro variations from sample to sample must be 
considered, rather than the actual values and uncertainty on R . Those variations come 
from possible thickness variations and are estimated to be smaller than 2%, in agreement 
with thickness measurements on Tl, T2, T3. Taking into account the uncertainty on the 
thermal resistance of the layers on SI, S2, S4 and the possible variation of the thermal 
resistance of the SiC>2 layers, figure ([3]) displays superlattice thermal conductivity variation, 
versus interface width. Clearly, SI and S2 exhibit close thermal conductivities whereas S4 
exhibits a significantly smaller thermal conductivity. 

Our data analysis neglects the boundary thermal resistances between the various lay- 
ers (line/insulating film, insulating film/superlattice, superlattice/substrate, insulating 
film/substrate). However, we may assume the boundary resistances are the same in the 
various samples. In that case, their contributions are partly cancelled out by the compari- 
son process between samples and test samples. Therefore, the thermal boundary effect may 
alter the absolute values we derived, but not the overall behavior of the thermal conductivity 
versus interface width. 
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FIG. 4: Aluminium concentration profile versus position along the growth axis z. (a) Perfect 

nn 

interface, (b) according t o llll . Il2fl. (c) do = 0.15 nm. (d) do = 0.45 nm. (c) and (d): dashed line: 
'erf profile according to [16| . Full line: average concentration in the monolayers. 



III. DISCUSSION 



Clearly, our experiments show that SI and S2 exhibit close thermal conductivities whereas 
S4 exhibits a significantly smaller thermal conductivity (30% decrease). 

In (GaAs)g(AlAs)5, the superlattice period d is much smaller than the thermal phonon 
mean free path I in the bulk constituents (d = 4 nm and I ~ 70 nm [9]). In this so called 
"short period superlattice" regime, we expect the interface defects to modify significantly the 
phonons mean free path and, as a consequence, the thermal conductivity. This is confirmed 
by molecular dynamics calculations. Direct comparison of our experimental results with 
quantitative models is difficult because existing calculations use simplified models for both 

nnn 

structure and interfaces |9J, HJJ, 112J • However, we may compare the transition layer at the 
interfaces involved in both models and experiments. Figure (jlj) is a sketch of the interface 
transition layers. In our experiments. SI. S2 and S4 arc characterized by an interface width 
do [lij]. From do, we may infer the foreign atoms concentration in the transition layer. In 
S2, do/x c = 0.5 (x c = 0.28 nm is the monolayer thickness) and the transition layer is two 
monolayers in thickness: the first and last monolayer at the interfaces. In both monolayers, 
the foreign atoms concentration is 15%. In S4, do/x c = 1.6 and the transition layer is 
four monolayers in thickness: the first and last two monolayers at the interfaces. In the 
monolayers the closest to the interfaces, the foreign atoms concentration is 33%. In the 
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next monolayers, the foreign atoms concentration is 10%. In existing molecular dynamics 
calculations, the transition layer is one monolayer in thickness. It is located at the last 
atomic layer in each superlattice layer and the foreign atoms concentration is 50%. A 60% 
reduction of the thermal conductivity is then inferred in (GaAs)7(AlAs)y [ill . [l^ |. On the 
basis of the above comparison, we conclude that our samples lie in a range where the interface 
damage should indeed have an impact on the thermal conductivity of our samples. 

Large changes of the total phonon mean free path and of the thermal conductivity have 
been predicted when diffuse scattering is introduced at the interfaces The largest sen- 
sitivity is achieved for small amounts of diffuse scattering. Therefore, the conductivity 
variation between SI and S2 is surprisingly small. The conductivity variation which is ob- 
served between SI and S4 is consistent with the order of magnitude predicted from molecular 
dynamics. However, the 30% variation which is measured appears to be quite small with 
respect the 60% variation predicted for a less severe interface damage. 

IV. CONCLUSION 

We have measured the cross-plane thermal conductivity of three short-period 
(GaAs)g(AlAs)5 superlattices, at room temperature. Due to different temperature growth, 
the samples exhibit different small scale roughness at the interfaces. The small scale rough- 
ness is characterized by an effective interface width do, on the basis of previous Raman 
scattering experiments. Within error bars, no variation of thermal conductivity is observed 
between d = and d = 0.15 nm, and a 30% decrease is observed between d = and 
do = 0.45 nm. On the basis of existing calculations, the sensitivity of the thermal conductiv- 
ity to the small scale interface roughness seems to be not as strong as expected. Additional 
experiments, as well as more realistic calculations must be performed to clarify this point. 
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